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Meissner screening masses in gluonic phase 
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A numerical analysis for the Meissner mass in the simplest gluonic phase (the minimal cylindrical 
gluonic phase II) is performed in the framework of the gauged Nambu-Jona-Lasinio model with cold 
two-flavor quark matter. We derive Meissner mass formulae without using the numerical second 
derivative. It is revealed that the gapless mode yields a characterized contribution to the Meissner 
mass. We also find that there are large and positive contributions from the tree gluon potential 
term to the transverse modes of gluons. It is shown that the simplest gluonic phase resolves the 
chromomagnetic instability in a rather wide region. 

PACS numbers: 12.38.-t, 11. 15. Ex, ll.30.Qc 



I. INTRODUCTION 

Quark matter at sufRciently high density and low tem- 
perature is expected to be in a color superconducting 
state driven by the BCS mechanism P, 0, Q. This is 
analogous to the electron Cooper paring in a supercon- 
ducting metal. However, quarks, unlike electrons, have 
color and flavor degrees of freedom as well as spin, so 
that the phase structure is quite rich. In nature, decon- 
fined quark matter might exist in the interior of neutron 
stars [4]. Thus the dynamics of the color superconduc- 
tivity has been intensively studied 

Bulk matter in compact stars should be in equilibrium 
under the weak interaction (/3-equilibrium), and be elec- 
trically and color neutral. The electric and color neu- 
trality conditions play a crucial role in the dynamics of 
the quark pairing 0, H, Hi- In addition, the strange 
quark mass cannot be neglected in moderately dense 
quark matter as in the compact stars. Then a mismatch 
6fi between the Fermi momenta of the pairing quarks is 
induced. 

As the mismatch 6ii increases, the conventional color 
superconducting state tends to be destroyed. Before the 
complete destruction, however, the Meissner mass of glu- 
ons turns to be imaginary in the gapped (2SC) and gap - 
less (g2SC) two-flavor color superconducting phases pi| : 
In the g2SC phase with the diquark gap A < Sfi 
the 8th gluon has an imaginary Meissner mass, while 
the Meissner masses for the 4- 7th gluons are imagi- 
nary also in the 2SC phase Sfi < A < ^/26fi. This 
chromomagnetic instability implies that there should ex- 
ist a more stable vacuum other than the 2SC/g2SC 
phase. Later a chromomagnetic instability was found 
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also in the three-flavor gapless color-flavor locked (gCFL) 
phase [111, [13, [l3| • One of the central issues in this field 
is to establish the genuine ground state for realistic val- 
ues of 6fi. Besides the gluonic phase [3, [lB| , a number 
of other candidates for the true vacuum have been pro- 
posed [li,[l3,[ll,[li,[2S[2ll[22,[2l[2i,[2i,[2^. 

Connected with the chromomagnetic instability, it was 
revealed that there appear tachyonic plasmons in the 4- 
7th and 8th gluonic channels [23|. It clearly shows that 
the physical vectorial excitations carry the instabilities 
and thus supports the scenario with gluon condensates 
(gluonic phase). 

It is also known that the physical diquark excitation 
(the diquark Higgs mode) suffers from the Sarma insta- 
bility 28j in the g2SC region, which corresponds to the 
negative mass squared of the diquark Higgs at zero mo- 
mentum. Furthermore, it was found that the diquark 
Higgs mode has a negative velocity squared < in the 
g2SC region [201 . A similar instability is also discussed 
in Refs. [30, [3l[ . This problem should be also resolved in 
the genuine ground state. 

In Refs. [3j[IB|j the Ginzburg-Landau (GL) approach 
in the hard dense loop (HDL) approximation was em- 
ployed in the vicinity of Sfj, ~ A/^/2. Outside the 
scaling region around (5/x ~ A/^/2, the self-consistent 
analysis by solving the gap equations and the neutral- 
ity conditions was recently performed in Ref. [s^: ft is 
shown that the gluonic phase is actually realized in a 
wide region of the parameter space and it is energeti- 
cally more favorable than the normal, 2SC/g2SC, and 
the single plane wave Larkin-Ovchinnikov-Fulde-Ferrell 
(LOFF) [11 [si, [ail phases. It is also found that the 
values of A and in the gluonic phase are significantly 
different from those in the 2SC/g2SC phase. It is notice- 
able that the values of the gluon condensate are large, 
say, C'(I00-250MeV) in the almost whole region where it 
exists. On the other hand, the values of the color chemi- 
cal potentials are relatively small. For the earlier works in 
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other approaches, see Refs. [35|,|36|. The extension to the 
model with nonzero temperature is studied in Ref. [s^ . 

In this paper, we examine whether or not the gluonic 
phase resolves the chromomagnetic instability. We derive 
formulae for the Meissner mass without using the numer- 
ical second derivative. It turns out that the gapless mode 
gives a special contribution to the Meissner mass. In the 
numerical analysis, we consider the gluonic phase with 
the simplest ansatz which is called the minimal cylin- 
drical gluonic phase II [l^, [s^. As a benchmark, we 
also analyze the single plane wave LOFF and 2SC/g2SC 
phases including the non-HDL corrections. 

We find that in the gluonic phase the tree gluon po- 
tential term yields positive and large contributions to the 
Meissner masses of the transverse modes of gluons. Ac- 
tually, in the minimal cylindrical gluonic phase II, the 
chromomagnetic instability is resolved in the weak and 
intermediate coupling region, 65.4MeV < Aq < 130MeV 
for = 400MeV and A = 653.3MeV, where ^ and 
A denote the quark chemical potential and the cutoff 
in the (gauged) Nambu-Jona-Lasinio (NJL) model, re- 
spectively. We here introduced the 2SC gap parameter 
Aq defined at (5/i = 0, which essentially corresponds to 
the four-diquark coupling constant in the (gauged) NJL 
model. In the intermediate and strong coupling region 
ISOMeV < Aq < 160MeV, however, the chromomag- 
netic instability occurs in the transverse modes of the 
4th and 5th gluons. Besides, in a small region around 
Aq — 150MeV, the squared Meissner mass for the trans- 
verse mode of the 8th gluon becomes negative. For the 
other modes, the chromomagnetic instability does not 
occur. On the other hand, the single plane wave LOFF 
phase resolves the chromomagnetic instability only in the 
region 64.9MeV < Aq < 80MeV. This is consistent 
with the results within the HDL approximation shown 
in Ref. js^]. Numerically, the 2SC phase suffers from 
the illness in 134.6MeV < Aq < 160MeV, whereas the 
g2SC phase does in Ao < 134.6MeV. We thus con- 
clude that the situation is definitely improved even in 
the simplest gluonic phase. If we consider more com- 
plicated gluonic phases such as the cylindrical gluonic 
phase I [ij, iT5| and /or the gluonic color-spin locked 
(GCSL) phase [1^, IS^], the chromomagnetic instability 
might be completely removed. 

The paper is organized as follows: In SecHIl the gauged 
NJL model is described. We also show the dynamical 
solutions of the gluonic, LOFF and 2SC/g2SC phases. 
In Sec lIII M we develop the formulae for the numerical 
calculation of the Meissner masses. We numerically an- 
alyze the Meissner masses for the gluonic, LOFF and 
2SC/g2SC phases in Sec HITBl Section HVl presents the 
summary and discussions. 



II. MODEL 



We study the gauged Nambu-Jona-Lasinio (NJL) 
model with two light quarks. We neglect the current 
quark masses and the ('0^)^-interaction channel. The 
Lagrangian density is given by 



C = V5(i|) -f /xo7°)^ + Ga 
1 



(V'^i£e"75V')(V'*ee"757/'^) 



(1) 



with 



(2) 

where s and e" are the totally antisymmetric tensors in 
the flavor and color spaces, respectively. We also intro- 
duced gluon fields AJ^, the QCD coupling constant g, the 
generators T° of SU (3), and the structure constants 
The quark field V' is a flavor doublet and color triplet. 
The charge-conjugate spinor is defined by ip'^ = C'ip'^ 
with C = Z7^7°. We do not introduce the photon field. 
On the other hand, the whole theory contains free elec- 
trons, although we do not show them explicitly in Eq. ([T|). 
In /3-equilibrium, the chemical potential matrix fiQ for up 
and down quarks is 



/XO = Ail - fleQc 



(3) 



with 1 = Ic 1/, and Qom = Ic (X" diag(2/3, — 1/3)/, 
where fi and are the quark and electron chemical po- 
tentials, respectively. (The baryon chemical potential /is 
is given by fiB = 3/i.) The subscripts c and / mean that 
the corresponding matrices act on the color and flavor 
spaces, respectively. Hereafter, we abbreviate the unit 
matrices, 1, Ic and 1/, if it is self-evident. By introduc- 
ing the diquark field ~ i-ip^ee^^^ijj^ we can rewrite 
the Lagrangian density ^ as 



£ = ^(z^ + ;xo7°)^-^-^$"[»#e"75V'''] 



(4) 



In the 2SC/g2SC phase, we can choose the anti-blue 
direction without loss of generality. 



($'■) = 0, ($5) = 0, = A, 



(5) 



where the diquark condensate A is real. In this basis, by 
imposing the color neutralit y c ondition, the color chem- 
ical potential /ig is induced |38j . We can interpret /ig as 
the vacuum expectation value (VEV) of the time compo- 
nent of the 8th gluon. 
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Let us define the Nambu-Gor'kov spinor, 



with 



* = 1 ,,c 



(6) 



The propagator inverse of in tlic 2SC/g2SC phase is 
given by 



with 



[G+r\p) = {po + ^l-5^iT3-^lslbh°-'f■p, (8) 



[Go]-\P) = ipo-^l + s^iT3 + ^lslbh''-1■P, (9) 



and 



A- = -zee^sA, A+ = -/''{A-y-f'' = -iee'^/^A, (10) 

where P'^ = {po,p} is the energy-momentum four vector. 
We also defined T3 = diag(l,— 1)/, 1;, = diag(0, 0, l)c, 
and 

The 2SC/g2SC phase is not the genuine ground state in 
the region Sii> A/\/2, because it suffers from the chro- 
momagnetic instability. A candidate to resolve the chro- 
momagnetic instability is the giuonic phase with gluon 
condensates 

Let us introduce the gluon condensates (A^) 7^ 0. 

When the space-component gluon condensates {A"-) ^ 
are incorporated into the theory, the time-component 
VEVs of the gluon fields other than the 8th one are gen- 
erally induced as well. We may interpret them as the 
color chemical potentials^ [soj : 

Ma = 9{A^), (a = 1, 2, • • • , 7), fis = ^5(^0) • (12) 

The propagator inverse of 'i' including the gluon con- 
densates is written as 



(13) 



[G+J-\P) ^ ipo + /xo)7" -1-P + 9{r)T\ (14) 
[G^J-\P) ^ {po - /xo)7" - 7 • P - girKT^'f- (15) 

In the fermion one-loop approximation, the bare effec- 
tive potential including both gluon and diquark conden- 
sates is given by 



V'barc 
eff — 



A^ 1 I r d^P 



AG A 4 



J^A /^Trln^-i 



(16) 

where we added the free electron contribution. Since the 
bare potential has a divergence, a counter term is re- 
quired. We take into account only differences of the free 
energies with and without the chemical potentials. We 
thus define the renormalized effective potential by 



T/i? T/barc 



Vet., 



(17) 



with the counter term, 



d'^P 
— — --Trln5'7"'" 
i(27r)4 9 



(18) 

In this prescription, even if we use the regularization 
scheme with the sharp cutoff A for the loop integral, we 
can remove artificial mass terms of gluons like A^A^. 

In general, we can reduce the 32 homogeneous gluon 
condensates to 25 ones [l^. However the general case is 
quite complicated and hence it is difficult to find the self- 
consistent solutions of the gap equations and the color 
neutrality conditions for the 25 VEVs at present. 

In this paper, we consider the minimal ansatz for the 
cylindrical giuonic phase II [H, , 



M8^^g(A«), B^g{Al). (19) 



In order to make the physical meaning of the gluon con- 
densates clearer, it is convenient to use the unitary gauge 
in which all gauge dependent degrees of freedom are re- 
moved. We shall fix the gauge as follows [31 : 



0, $9 = 0, Im$'' = 0, 



(20) 



and 



^ In an appropriate basis, the color chemical potentials can be 
reduced only into fi^ and , because the color chemical potential 
matrix fic = MaT" (a = 1, 2, ■ • ■ , 8) is hermite and traceless. In 
this case, however, the basis for the diquark field changes from 
Eq. (0. 



A! = 0, = 0, Al = 0. 



(21) 



As a benchmark, we also consider the single plane- wave 
2SC-L0FF phase [ii,!!!!!!. 



($'■) = ($9) ^ 0, ($'') = Ae-2*«-^ , 



(22) 



60 80 100 120 140 160 
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FIG. 1: The dynamical solutions for the minimal cylindrical 
gluonic phase II. The bold solid, bold dashed, thin solid, thin 
dashed and thin dot-dashed curves represent the values of B, 
A, 5fi, ^3 and /is, respectively. The values A — 653.3 MeV 
and fj, — 400 MeV were used. We also took as = 1. 

Introducing the quark field "0' = e^'^'^Tp with the color- 
singlet phase, we can erase the x-dependent phase of the 
LOFF order parameter and instead, the x-independent 
color-singiet term -tpyqip is induced in the kinetic term for 
quarks. Although the vector (f is also gauge equivalent to 
the condensate {Ag), there is subtlety with respect to the 
tree gluon kinetic term: Notice that the tree gluon po- 
tential does not give any contribution to the free energy, 
while it is relevant to the Meissner masses for the trans- 
verse modes of the 4- 7th gluons, because does not 
commute to T^~'^ and thereby the corresponding Meiss- 
ner masses are of the order of q{= \q\). The point is that 
the values of q are large in a wide range of the parameter 
region where the LOFF phase exists, as we will see below. 
In order to avoid confusion, we may use the color-singlet 
transformation or we simply do not incorporate the tree 
gluon potential term into the LOFF phase in any case. 

The dynamics of the minimal cylindrical gluonic 
phase II is analyzed in Rcf. 32] by solving the gap equa- 
tions and the neutrality conditions in a self-consistent 
way. We explicitly show the results for the gluonic, sin- 
gle plane wave LOFF and 2SC/g2SC phases in Figs[T]- 
[31 respectively. In the analysis, we took realistic values 
fi = 400MeV and A = 653.3MeV. We also converted 
the four-diquark coupling constant Ga to the 2SC gap 
parameter Aq defined aX S/i = and varied the values 
of Aq from the weak coupling regime (Aq ^ 60 MeV) to 
the strong coupling one (Ao ~ 200 MeV). For the gluonic 
phase, it is required to specify the value of as[= (?^/(47r)], 
although the results for the minimal cylindrical gluonic 
phase II are not sensitive to the choice of the values of 
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FIG. 2: The dynamical solutions for the single plane wave 
LOFF phase. The bold solid, bold dashed, thin solid and 
thin dot-dashed curves represent the values of q (= \q\), A, 
Sjj, and ^8, respectively. The values A — 653.3 MeV and 
fj, = 400 MeV were used. 

as ;32]. We here took — 1. 

While the neutral normal phase without the diquark 
condensate always exists, the neutral gluonic, LOFF, 
g2SC and 2SC phases do only in the regions, 

65.4MeV < Ao < 160McV, (gluonic) (23) 
64.9MeV < Aq < 138MeV, (LOFF) (24) 
92.2MeV < Ao < 134.6MeV, (g2SC) (25) 

and 

Ao > 134.6MeV, (2SC) (26) 

respectively. 

The analysis for the free energies has been done in 
Ref. (3^ : The normal phase is realized in the weak cou- 
pling regime with Ao < 64.9MeV. While the single 
plane wave LOFF phase is energetically most favorable 
only in the narrow region 64.9MeV < Ao < 67MeV, 
the minimal cylindrical gluonic phase II is stabler than 
the LOFF and 2SC/g2SC phases in the wide region 
67MeV < Ao < 160MeV. In the strong coupling regime 
with Aq > 160MeV, the 2SC phase is reahzed. 

For the numerical calculation of the Meissner masses, 
we use the solutions shown in Figs[T]-[31 It is noticeable 
that the condensate B is large in the almost whole region 
where the gluonic phase exists. (See Fig. [TJ) This feature 
is crucial for the Meissner masses in the gluonic phase, 
as we will see in the next section. 
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where we defined the tree contribution 



_|_^2 jaiaa2 jaiba^ ^^v j^a2 j^a^ ) 



and the vertex 



1 ds-^ 

p^a = g - 



dAl 



7^T'^ 
_^A'(T'^)r 



(28) 



(29) 



In Eq. (P7| . "c.t." denotes the counter term and we used 



FIG. 3: The dynamical solutions for the 2SC/g2SC phase. 
The solid, dashed and dot-dashed curves represent the values 
of A, and /xg, respectively. At Ao = 134.6 MeV, the g2SC 
phase turns into the 2SC one. The values A = 653.3 MeV 
and [1, — 400 MeV were used. 



III. MEISSNER MASSES IN GLUONIC PHASE 
A. Formulae 

Let us derive formulae for the numerical calculation of 
the Meissner screening mass. 

The squared Meissner mass can be expressed through 
the second derivative of the effective potential: 



dA-dA<l 



' 2 



5^ r d^p 
~2 J i{2^Y 



Tr 



i{2T:y 



-Tr 



SgV^SgT 



vb 



^-^iSgS-') 



^"g c— 1 I c g 



(30) 

for X =^ j4^, A,(5/i, and linearity of 5*^"^ with respect to 



all variables, i.e.. 



0. We also abbreviated the 



•' dXdY 

bracket (• • •) for the gluon condensates. 

For the numerical calculation, it is useful to diagonalize 
Sg^^ and/or Sg. Although the propagator inverse in 
Eq. (fT3|) is a 48 X 48 matrix in the flavor, color, spinor and 
Nambu-Gor'kov spaces, we can block-diagonalize 5*^"^ in 
the flavor and chirality spaces. This technique reduces 
our labour. 

Let us transform the propagator inverse as follows; 



Y> 

i£75 



, (27) with 



I , (31) 



S7\P)=pol + Hg, Hg = -S^IT3 



-e''A 



e'-A 

-p. - (1^ + 7°7 • p- 7°7 • A^ 



(32) 



r 



where 



y, /^c 



gA^T", A = gA^T" 



(33) 



flavor-diagonal. Since the inverse of Eq. ([31]) yields the 
expression for the propagator. 



We here decomposed Sg ^ into the diagonal po-part and 
the "Hamiltonian" Hg. (One can check easily hermitic- 
ity of Hg.) Notice that the flavor dependence of 
exists only in the first term of Hg and therefore is 



7" 
ie-f5 



(34) 
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the second derivative then reads 

.2 , ^4p 



dA-dAl 



tree 



' 2 



9_ 

2 j i{2^Y 



Tr 



i(27r)' 



-Tr 



SgV^SgV 



vb 



, (35) 



One can easily show that the trace over the spinor space 
in Eq. ([BS]) is the sum of the two parts. 

In virtue of hermiticity, we can diagonahze Hg and Sg 
by a unitary matrix [/, 

Hg = UH^,^^U\ iJdiag - diag(i?[, El, • • • , El), (38) 



with 



7 







«£75 



1 







ie^°^b 



-7M^0(r'^)T 



(36) 



(37) 



Since the current quark masses are ignored, the theory is 
chiral invariant. Thus we can decompose the vertex and 
also the propagator into the right and left-handed parts. 



Sg — [/diag 



1 



1 



Po + El'po + Ef 'po + E-^ 

(39) 

where r = ± for — ±1 and EI2 ... „ denote the en- 
ergy eigenvalues. It is not difficult to find numerically 
the energy eigenvalues and the unitary matrix by using a 
standard method. Noting that the integrand of Eq. ({551) 
contains the po-dependence only in Sg with the expres- 
sion (|39p , we can explicitly perform the integral over pq 
and thereby obtain 



dA-dAi 



tree 



2 ^ ^ 7 (27r)3 EJ - E 



d^p ^(i^r)-^(^(^tfM<^c/).,(t/tf''''c/), 



' 2 



E E 

T=± et=e: 



d^p 
(2^ 



6iEl){U''r'"'U)ij{U''r''''U)ji - (counter term) . 



(40) 



We can derive similar formulae for other second deriva- 
tives, 



dMAl ' 



etc.. 



(41) 



It is also straightforward to extend the formulae to the 
version with a finite temperature. 

The formula l|40p has an advantage over the numerical 
derivative of the effective potential: The numerical sec- 
ond derivative requires a quite precise calculation for the 
free energy and thereby it takes a long time. On the other 
hand, we can reach a sufficiently accurate result via (HDl) 
in a reasonable time. Furthermore, in the expression of 
PO)l , it is clear that the contributions of the gapped and 
gapless modes to the Meissner masses are quite differ- 
ent. (Compare the second and third terms in Eq. ([M)) .) 
This might help us to understand why the existence of 
the gapless modes yields a sudden change of the squared 
Meissner mass for the 8th gluon at the border of the 2SC 
and g2SC phases p^ . 



B. Numerical analysis 

Before the numerical calculation, we describe several 
features of the Meissner masses in the minimal cylindrical 
gluonic phase II. 

In this phase, the symmetry breaking structure is [l5l | 



[5f/(3)c]local X [/(l)e,n X 50(3),ot 

X 50(2)rot, 



^^(1) 



(42) 



where the unbroken ?7(l)om is connected with the new 

electric charge Qcm = Qcm ~ ^-^^ ^ -^"^^ ^^"^ rotational 
symmetry breaking leads to different Meissner masses for 
the transverse (/i, v = x,y) and longitudinal (/i, v = z) 
modes. It is thus convenient to define 



[MYab 



(MYab 



dA%dAl 
dA'idA'l 



dA^dAl ' 



(43) 
(44) 



Since we took the unitary gauge (|21|). the squared Meiss- 
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FIG. 4: The squared Meissner masses of the transverse mode 
of the 4th gluon in the unit of Mg [= 4as/i^/(37r)]. The 
bold sohd, dashed and dot-dashed curves are for the minimal 
cylindrical gluonic phase II, the single plane-wave LOFF and 
the 2SC/g2SC phases, respectively. The values A = 653.3 
MeV, fj. = 400 MeV and = 1 were used. 

ner masses for the physical degrees of freedom are 

(M^)^!^ = (M^)^2^ = (M^)^5, (45) 

and 

iM')lf, (M2)^•^ (M^)^,, iM^)l-\ (46) 

where the relations (j45p hold owing to the unbroken 

C/(l)om symmetry. For the transverse modes of the 3rd 
and 8th giuons, it turns out that there exists a large mix- 
ing term (M^)|g, so that we define the diagonal mass- 
squared terms for them, 

(^^')^3,diag, (^')^8,diag- (47) 

Do there exist two Nambu-Goldstone (NG) bosons 
connected with the symmetry breaking SO (3) rot ^ 
50(2)1.01? This is nontrivial because Goldstone's theo- 
rem for relativistically invariant theories is not necessar- 
ily valid in noninvariant systems [iO, |4l| . We find that 
the answer is formally "yes" in this case, as we will see 
below. 

The point is that the rotational symmetry is sponta- 
neously broken only by the condensate (A^) ^ in the 
minimal cylindrical gluonic phase II. Therefore, before 
taking the z-direction, the effective potential should de- 
pend on the S0{3)rot invariant 

(^')'' (48) 

i—x^y^z 



FIG. 5: The squared Meissner masses of the transverse mode 
of the 8th gluon in the unit of Mg [= 4asfJ,^ /{3tt)]. The bold 
solid and dot-dashed curves are for the minimal cylindrical 
gluonic phase II and the 2SC/g2SC phases, respectively. The 
values A — 653.3MeV, n — 400MeV and cts = 1 were used. 
In the g2SC phase, the squared Meissner mass monotonously 
decreases and goes to minus infinity at Ao = 134.6MeV, al- 
though it is not explicitly shown here. 

i.e., 

Kf = Kf(S). (49) 

We then find 

Since we take the direction {A^ y) = 0, (A^) ^ and the 

gap equation for (A^) yields = 0, we formally find 
that the squared Meissner mass for the transverse mode 
of the 6th gluon is vanishing, 

(M')i; = 0. (51) 

It implies that A^ y correspond to the two NG bosons. 

A crucial difference between the gluonic and 2SC/g2SC 
phases is the existence of the tree gluon potential term. 
Although the effect is negligible for the free energy in the 
minimal cylindrical gluonic phase II, it is quite important 
for the Meissner masses. Neglecting the suppressed terms 
~ O {fj,^) , O {113 ^s) , (^(Ms)' ^e obtain the tree terms of the 
squared Meissner masses: 

(M2)n,22 ^ (M')r4,55 ^ iM^)l3 = (52a) 
{An^,^B^ (A^^)i, = Js^ (52b) 

{M')Is--^B\ (52c) 
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FIG. 6: The squared Meissner masses of the longitudinal 
mode of the 6th gluon in the unit of Mg [= 4Qs/i^/(37r)]. 
The bold solid, dashed and dot-dashed curves are for the 
minimal cylindrical gluonic phase II, the single plane- wave 
LOFF and the 2SC/g2SC phases, respectively. The values 
A = 653.3MeV, fj, = 400MeV and a, = 1 were used. 

Thus the transverse modes except for (M^)|"g and 
diag have the positive and large contributions of 
the order of B^. (For the values of B, see Fig. [T]) This 
is one of the reasons why the Meissner masses tend to be 
positive compared with those in the 2SC/g2SC phase. 

We also note some features of the Meissner masses for 
the single-plane wave 2SC-L0FF phase. The Ajj~^ glu- 
ons should be massless and the relations 

(M^)^!^ = (M2)^f = (M^)^^^ = iM^)^f (53) 

hold because of the unbroken SU{2)c gauge symmetry. 
In addition, similarly to (I5ip . we formally obtain 

{M')Js^O- (54) 

Let us now turn to the numerical analysis of the Meiss- 
ner masses. 

We depict the results in FigsHHSlin the unit of Mg[= 
4Q;sM^/(37r)]. In the analysis, we used — 400MeV, A = 
653.3MeV and = 1. 

For the minimal cylindrical gluonic phase II, the 
squared Meissner masses (M2)f4 = (M^)!^ are positive 
in the region 65.4MeV < Ao < 130MeV, while it suffers 
from the chromomagnetic instability in 130MeV < Aq < 
leOMeV. (See Figgl) We also find that (M2)|g becomes 
negative in the small region around Aq ^ 150McV. (See 
FigEl) For the other modes, however, the chromomag- 
netic instability does not occur as shown in FigsUHHl (Af- 
ter the diagonalization of (M'^)'^^, {M'^)38 and (M2)^g, 
the instability in (Af^)|"g is converted into {M'^)^^ diaz^ 



FIG. 7; The squared Meissner masses of the longitudinal 
mode of the 8th gluon in the unit of Mg [= 4as/i^/(37r)]. 
The bold solid, dashed and dot-dashed curves are for the 
minimal cylindrical gluonic phase II, the single plane- wave 
LOFF and the 2SC/g2SC phases, respectively. The values 
A = 653.3MeV, = 400MeV and as = 1 were used. As in 
FiglU the curve for the g2SC phase below the figure is cut ofi'. 



because we define the diagonalized squared masses as 
(^')i3,diag < (^')?8,diag-) ^ q^ite noticeablc that 
there exist spikes and valleys around Aq ^ lOOMeV and 
Ao - 150MeV in FigsUH 

How about the sensitivity of the Meissner masses on 
Q!s? Although the dynamical solutions of A, B, (5/x, 
/i3 and /is in the minimal cylindrical gluonic phase II 
are almost independent of [S^l , the Meissner masses 
for the transverse modes can be sensitive. Note that 
the one-loop contributions are proportional to as and 
thus the influence of the tree contributions ([5^ is rel- 
atively stronger (weaker) as the values of ag decrease 
(increase). For example, (M^)44 becomes negative at 
Ao = 140, 130, 120MeV for = 0.85, 1.0, 1.15, respec- 
tively. For the other transverse modes in Eq. (|52p . there 
should appear similar sensitivities. On the other hand, 
for the longitudinal modes, the ratio (M^)^j/M^ is in- 
sensitive to as, because the tree contributions are sup- 
pressed. 

For the single plane wave LOFF phase, like in the min- 
imal cylindrical gluonic phase II, the transverse modes of 
the 4- 7th gluons are the most problematic. However the 
chromomagnetic instability occurs in the earlier region, 
80MeV < Ao < 138MeV. (See FigH) All the other 
modes also suffer from the chromomagnetic instability in 
the end of the LOFF phase around Ao 130MeV. (See 
FiglS] and [71) This is different from the situation in the 
minimal cylindrical gluonic phase II. We also note that 
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FIG. 8: The squared Meissner masses for several gluonic 
modes in the minimal cylindrical gluonic phase II in the 
unit of [= 4as/iV(37r)]. The values A = 653.3MeV, 
fi — 400MeV and a, = 1 were used. In the right bottom 
figure, we showed the diagonalized mass squared terms for 
the transverse modes of the 3rd and 8th gluons. 



these results are consistent with the analysis based on 
the HDL approximation [33l[33|. 

The 2SC/g2SC phase has the chromomagnetic in- 
stabihty numerically in the region 92.2MeV < Aq < 
IGOMeV. The results agree with those in Appendix B 
in the second paper of Ref. [lo| (, see also Refs. [H,!!!]). 

In conclusion, although the minimal cylindrical gluonic 
phase II does not completely remove the chromomagnetic 
instability, it cures the situation in a wide region. 

The comments about the massless modes are in or- 
der. The numerical calculations including the non-HDL 
effects do not necessarily reproduce the vanishing Meiss- 
ner masses. This fact is known even in the 2SC phase, 
i.e., the Aj^~^ gluons acquire the non-HDL contribu- 
tions like log A^/A^ in the sharp-cutoff regularization 
scheme [H, . In order to settle this problem, a more 
sophisticated regularization scheme is required. It will 
be studied elsewhere. 



IV. SUMMARY AND DISCUSSIONS 

We analyzed the Meissner screening masses in the 
simplest gluonic phase (the minimal cylindrical glu- 
onic phase II) as well as the single plane wave LOFF 
and 2SC/g2SC phases. We derived the formulae for 
the Meissner masses without any help of the numerical 
derivative. It was found that in the formulae the gap- 
less mode makes the contribution characterized by the 
Dirac's (5-function. We showed that the simplest gluonic 



phase removes the chromomagnetic instability in the re- 
gion 65.4MeV < Aq < ISOMeV, whereas the single plane 
wave LOFF one does in 64.9MeV < Aq < 80MeV. We 
here took A = 653.3MeV, fj. = 400MeV and = 1. The 
2SC phase does not have the chromomagnetic instability 
in the strong coupling regime Aq > IGOMeV. Incorpo- 
rating the analysis of the free energy [32|, we conclude 
that in the region 67MeV < Aq < ISOMeV the simplest 
ansatz for the gluonic phases works more nicely than the 
single plane wave LOFF and g2SC phases. On the other 
hand, the single plane wave LOFF phase is energetically 
more favorable and also resolves the chromomagnetic in- 
stability only in the window 64.9MeV < Aq < 67MeV. 



Furthermore, we found the noticeable behaviours of 
the squared Meissner masses in the minimal cylindri- 
cal gluonic phase II around Aq ~ lOOMeV and Aq ^ 
150MeV. (See FigsgHH) The 2SC/g2SC phase also has 
the similar behaviours: Notice that there appears the 
abrupt change of the Meissner mass for the 8th gluon 
from the g2SC side to the 2SC one and that the values of 
the squared Meissner mass for the 4- 7th gluons have the 
valley at the phase transition point from the g2SC phase 
to the 2SC one, as shown in FigsHHTl These similarities 
might suggest that new gapless modes in the minimal 
cylindrical gluonic phase II appear around Aq ~ lOOMeV 
and Aq ~ ISOMeV. We also note that the value of B 
takes its maximum around Ag ~ lOOMeV and that the 
relation i? ~ is satisfied around Aq ~ ISOMeV. (See 
FigUI) These facts might be significant. The dispersion 
relation for quarks in the gluonic phase will be performed 
elsewhere. 



There still exists a chromomagnetically unstable re- 
gion. However it should be noticed that we examined 
only the simplest ansatze for the LOFF and gluonic 
phases in this paper. The multiple plane wave LOFF 
phase may completely remove the chromomagnetic in- 
stability in the whole parameter region [ly]. More 
involved gluonic phases can also resolve the instability: 
Since there appears the illness in the transverse mode of 
the 4th gluon, the GCSL phase with the gluon conden- 
sates Us = \/3/25(Ag) and K = g{Al) = g{Al) [H,!!! is 
hopeful, for example. An important point is that the free 
energy for the GCSL phase is slightly lower than that for 
the minimal cylindrical gluonic phase II Ii32il . Inhomoge- 
neous gluonic phases are also interesting [23, . 



Independently of the chromomagnetic instability, the 
Sarma instability for the diquark Higgs mode should be 
removed as well. This problem will be considered else- 
where. 
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